Structure formation on the brane: A mimicry 
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We show how braneworld cosmology with bulk matter can explain structure formation. In this 
scenario, the nonlocal corrections to the Friedmann equations supply a Weyl fluid that can dominate 
over matter at late times due to the energy exchange between the brane and the bulk. We demon- 
strate that the presence of the Weyl fluid radically changes the perturbation equations, which can 
take care of the fluctuations required to account for the large amount of inhomogeneities observed in 
the local universe. Further, we show how this Weyl fluid can mimic dark matter. We also investigate 
the bulk geometry responsible for the scenario. 

PACS numbers: 04.50.+h, 98.80.-k, 98.80.Cq 



I. INTRODUCTION 

£N| ■ Observations suggest that the expanding universe is homogeneous and isotropic at scales larger than 150ft. -1 Mpc. 
But the large amount of inhomogeneities observed in the local universe needs sufficient and convincing explanation. 
The usual cosmological models based on the standard Friedmann equations require that the baryonic matter fluc- 
tuation 6 > 1 today, implying S > 10~ 3 at the time of recombination. This is in direct contradiction with cosmic 
microwave background observations by over an order of magnitude. This limitation of standard cosmology leads to 
an inevitable prediction of dark matter p|. CMB Observations @ suggest that f^baryon^ 2 ~ 0.02 i.e., only about 
4% of cosmic density is baryonic. Hence dark matter, if it exists, has to be the dominant nonbaryonic component, 
contributing to as much as l/3rd of cosmic density It is proposed that its nonbaryonic nature helps it decouple 
' from radiation, resulting in a growth of structure that starts much before the hydrogen recombination. But there are 
q-i several problems associated with dark matter, the most pronouncing of which being its ill-response to detection by 
series of experiments. Though the experiments of DAMA group Q have reported in favor of its existence in form of 
(3J[v weakly interacting massive particles (WIMP) [Bj], similar searches by a number of other groups, such as CDMS [||, 
CRESST 7], EDELWEISS y| have yielded negative results (for recent results, see Q)- So, there arise questions such 
. £h [ as whether dark matter exists at all or it is the gravity sector, rather than the matter sector, that needs modifications. 

This leads people to consider modified gravity theories, eg, modified Newtonian dynamics (MOND) ^(Jj bifurcating 
^ ' gravity phantom cosmology etc. 

One such modified gravity theory is the braneworld gravity 01 which has opened up new avenues of explaining 
the observations with the help of a modified version of the standard Einstein equation [TJJ. In this scenario, when the 
bulk consists of matter, the bulk metric for which the FRW geometry on the brane is recovered, is given by a higher 
dimensional generalization of the radiative Vaidya black hole |15| that exchanges energy with the brane |l6t Il7| . In 
presence of bulk matter the modified Einstein equation on the brane reads [l6| 

k 2 

= -Ag^u + k 2 T^ v + 6— - £ M „ + (1.1) 

where S^ v , £ M „ and T^v are the quadratic contribution from brane energy-momentum tensor, the projected bulk Weyl 
tensor and the projected bulk energy-momentum tensor on the brane respectively. 

Incorporating all the braneworld corrections, one can conveniently express the Friedmann equations on the brane 

as 

H 2 = ^p cff +^-4f (1-2) 
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// = -_(p<= ff + p cff ) + 4 r (1.3) 



2 xr r ' a 
where the effective density and pressure are given by |l8j | 

f. f+ ' + W (1.4) 
2A a* 

The term p* = C(t)/a A is the combined effect of and .F^, and is called the "Weyl fluid", that supplies an 
additional perfect fluidlike effect to the usual brane fluid. Further, the Weyl parameter C(t) is related to the black 
hole mass, that results in energy-exchange between the bulk and the brane so that the individual matter-conservation 
is no longer valid but the total mass-energy of the bulk-brane system is now conserved. Consequently, the matter 
conservation equation on the brane is no longer a sacrosanct equation. Rather, it is modified to the nonconservation 
equation 

p + 3-(p + P ) = -2iP (1.6) 
a 

where ip = dm/dv, with m(v) the mass of the Vaidya black hole. It shows how the brane either loses (ijj > 0) or 
gains (jp < 0) energy in exchange with the bulk black hole. To a braneworld observer, the term ip is the quantitative 
estimate of the brane-projection of the bulk energy density. 

The goal of the present article is to show that the bulk-brane energy-exchange results in the growth of the Weyl 
fluid at late times. Following Newtonian analysis of perturbations from gravitational instability, which is the simplest 
yet logical analysis of gravitational perturbations, we demonstrate that the Weyl fluid can mimic dark matter in 
explaining structure formation. That the braneworld gravity can be a very good alternative to dark matter in 
astrophysical contexts of clusters and galaxies was proposed in [2(1 \2l\ . In this article we address the cosmological 
sector of dark matter. 

The plan of the paper is as follows. In Section-II, we obtain the perturbation equations for the effective perfect 
fluid on the brane. Section- III is devoted to the solution of the effective perturbation equations with the help of Weyl 
fluid that mimics dark matter, followed by a comparative analysis with the other dark matter models from standard 
cosmology as well as modified gravity theories. We construct the bulk geometry for this setup in Section-IV. Finally, 
we summarize our results and discuss some open issues. 



II. EFFECTIVE PERTURBATION EQUATIONS 

Since our focus is on the late time behaviour in the matter- dominated era, we restrict ourselves to the analysis 
of the zero brane cosmological constant scenario. Hence, the equations of hydrodynamics that involve the quadratic 
brane correction and the Weyl fluid correction to the brane perfect fluid, are 



dp 



,ctr 



Ot 



,cff 



>off ) = 


(2.1) 


- ^$ cff 


(2.2) 


47rG j o off 


(2.3) 



V 2 $ eff 

where lf cS is the velocity field in the effective perfect fluid. It should be noted that the term $ cff is not the usual 
Newtonian potential but the effective gravitational potential which is the resultant effect of the Newtonian as well as 
the relativistic potential. The later plays a crucial role in the braneworld context and has been discussed in details in 
pcj. We next consider small perturbations 

p cS (lt, t) = (f s (t)(l + 5 cS (^, t)) (2.4) 
$ eff (-^, t) = <$>f + (jf s (2.5) 

where p eS {t) and $q S are respectively the unperturbed effective density and effective potential and 5 eS and <p eS are 
their corresponding fluctuations. It is worthwhile to mention the significant difference of the density fluctuation of 
the braneworld cosmology from that of the standard cosmology. In the standard cosmology, 5 is the fluctuation of 
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baryonic matter only. On contrary, in the braneworld cosmology, <5 eff if the sumtotal of the fluctuations of baryonic 
matter and of the contribution from braneworld corrections. 

We proceed by expressing in terms of comoving coordinates and neglecting terms of second or higher order, which 
results in the following set of simplified perturbation equations 



ell' 



at 



—►off n 
U =0 



dt 



! + ^e ff = _l^_l^e fi 
a a p ctt a 

V^ cff = 4nGa 2 ff s S cS 



and express the solution in terms of Fourier transform 



(2.6) 

(2.7) 
(2.8) 



5° s (t, i) = £*f(*) 



ST(t) = y I <5 C «(^, t) e 



a k .it 



d 3 l? 



(2.9) 
(2.10) 



Further, assuming the effective pressure to be a function of the effective density alone, the equations of hydrody- 
namics now transform into a linear perturbation equation for 8^ : 



where c 2 eS = fa 



p+p 



d 2 Sf 
dt 2 



4p* 



2 ^ d Jt-{^Gf s -{cl ti k/af)5f = 
l 



p+A "T" 9(p+p)(l+p/A) 



a dt 
1 H 



(2-11) 



4p* 

3(p+p)(l+p/A) 



is the effective sound speed squared H 



III. SOLUTIONS WITH THE WEYL FLUID 



The above perturbative analysis can account for the required amount of gravitational instability if the Weyl density 
redshifts more slowly than baryonic matter density, so that even if it starts from a small initial value, it can eventually 
dominate over matter. Now, the nature and evolution of the Weyl fluid is governed by Eq l|1.6l) via the 4D Bianchi 
identity WG^ = 0. This gives 

p*+A-p* = Q (3.1) 
a 

where Q is a coupling term which can be calculated, on principle, if the projected bulk energy density ip is known. 
But in practice, since no one can fix the exact bulk geometry a priori, one has to take an ansatz for Q so far as it 
is physically reasonable and consistent with the brane equations. One such ansatz has been considered in 22] for a 
dilaton field in the bulk. Let us here take an ansatz Q = aHp* (a > 0), for which the Weyl fluid behaves like 

p* oc —r- — r (3.2) 

Thus the Weyl parameter is given by C(t) = Co a a (t), where Co is its initial value at the matter-dominated epoch. 
Obviously, the Weyl fluid is strictly radiationlike only if a = 0, i.e. for matter-free bulk scenario. But for the bulk 
with matter, the nature of the Weyl fluid depends on the coupling strength a. The more the coupling strength a 
(within the range 1 < a < 4), the more the dominance of the Weyl fluid over matter. Hence the Weyl fluid can mimic 
dark matter for a f=a 1. However, calculations of CMB anisotropies from the present model may lead to a better 
quantitative estimation for a. From Eq l|1.6|) . a > =>■ -0 > reveals that the brane loses energy to the bulk black 
hole. The increase of black hole mass is felt by a braneworld observer through the projected bulk energy density. 
Consequently, it results in the growth of the Weyl density at the expense of brane energy. 

We are now in a position of dealing with the perturbation equation (|2.11l) . This involves the fluctuation of p eS 
which is a sumtotal of three quantities given by Eq (|1.4|) . Of them, the quadratic correction term p 2 /2A comes into 
play at the physics of early universe such as during inflation (where p> A) [23j] but it contributes very little at the 
present era since p -C A > (lOOGeF) 4 . Hence, for all practical purpose, the effective density at late times can be 
approximated as 



P cn *p + p* (3.3) 
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What turns out from the above equation is that along with the usual matter density, here we have an additional 
(Weyl) density contributing to the total density that governs the perturbation equation l|2.11[l . Separating the baryonic 
(matter) part from the nonbaryonic (Weyl) part of Eq l|2.11l) now yields two wave equations 

= 4nGp B S B + 4mGp*5* (3.4) 
= iirGp*5* + AirGp B d B (3.5) 



d 2 6 B 


a dSn 


dt 2 


a dt 


d 2 5* 


a dd* 


dt 2 


a dt 



where 5b and 5* are the fluctuations of baryonic matter and Weyl fluid respectively and we have neglected the term 
involving sound speed because of the growing fluctuations. Consequently, with £lb <C f2*, the relevant growing mode 
solution for Eq l|3.5|l . as a function of the redshift, turns out to be 

S*(z) =6*(Q)(l + z)- 1 (3.6) 

Substituting the above expression in the fluctuation equation (|3.4|l of baryonic density gives 

^ + 2^ = 4 7 rGp^(0)(l + z)- 1 (3.7) 

Further, we know that though the universe evolves differently at early times, the standard cosmological solution 
for the scale factor is recovered in RS-II type brane- world gravity at late times [la. l24j . So, the late time behaviour 
for a spatially flat brane is given by 

g N 2/3(iu+l) 



a(t) = ^-H t) (3.8) 

With this scale factor and considering f2* i=s 1 at present time, Eq (|3.7J) now takes the form 

a 3/2 d ( a -i/ 3 dM + 2 dSs = 3 
da \ da J da 2 

A typical solution for the above equation is given by 

5 B (z) = 5*(z)(l-^±f) (3.10) 
V 1 + z N J 

where we have used the standard relation of the scale factor with the redshift function a oc (1 + z)" 1 . Note that 
though the above relation is a look-alike of the standard cosmological relation, physically it is completely different. 
Unlike the usual dark matter fluctuation, here 6* is the fluctuation of Weyl density that arises naturally in braneworld 
context. 

Let us now analyze some of the basic features of the present model. Eq l|3.10|l reveals that at z — » z/v, the baryonic 
fluctuation 5b — ► while 5* remains finite. This implies that even if the baryonic fluctuation is very small at a redshift 
of zn ~ 1000, as confirmed by CMB data @, the fluctuations of the Weyl fluid had a finite amplitude during that 
time. At z zn the baryonic matter fluctuations are of equal amplitude as the Weyl fluid fluctuations. This explains 
the structures we see today. Further, in this perturbative analysis, no extra matter (e.g. dark matter) has to be put 
by hand in order to explain the structures we see today. 

The effective equation of state parameter is given by 



..off _ = P + p(p + 2p)/2X + C(t)/ 3a i 
p eff p + p 2 /2X + C(t)/a 4 



off P h* ~ r\l ' III v ~ ^ V // foil\ 
W = Toff = „, J/OV , ^W~4 ( 3 " U ) 



which, in the matter-dominated era, can be approximated as 

wCff * 3(l + ki-«) (3 - 12) 

Clearly, it bears significant difference from the equation of state of cold dark matter (w = 0). In this sense the term 
'mimicry' may sound a bit misleading, though here we mean that it is the perturbative nature of evolution that is being 
mimicked by the Weyl fluid. Further, w cS does not cross the phantom divider line (w < —1) [2^| at least during the 
matter-dominated era we are interested about. Hence this model is favored by SNLS 26]. Moreover, the braneworld 
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scenario provides us with a new window for the cosmic coincidence problem. Here the evolution of the universe may 
not be a cosmological constant (or dark energy) effect at all, rather an outcome of the leakage of gravitational signal 
into the extra dimension. Hence the focus now shifts from the coincidence problem to the question: how can the 
late-acceleration of the universe be explained by the modified Friedmann equations ? 

A comparative study of our model with the other dark matter models and modified gravity theories is instructive. 
Eq (|3.12|l reveals that our model is distinct from the widely popular ACDM model p| so far as the equation of state 
is concerned. However, quite surprisingly, the scale factor bears strong similarity with that of the matter-dominated 
phase of ACDM. Hence, though they have some common features, whether or not the present model will eventually 
evolve into ACDM cosmology at very late times needs further investigation. As of now, this is an open question 
as obvious from the preceeding discussions. The alternative gravity theory modified Newtonian dynamics (MOND) 
[Tof . is based on modifying Newton's law. Though Bekenstein has recently proposed a relativistic MOND p^|. a 
more convincing and satisfactory relativistic version bearing important features such as lensing is yet to come. In 
comparison, the braneworld model is based on purely relativistic idea and hence is potentially more advanced. Another 
alternative cosmology is the bifurcating theory [llj where the effective Lagrangian bifurcates into several branches, 
thereby giving a possibility of unifying dark matter and dark energy. However, here a scalar field is still required to 
model dark matter. The model discussed in the present article do not need any such scalar field. Here the bulk-brane 
geometry plays the trick. Phantom cosmology |l2( can also provide an unified description but it has a wrong sign in 
the kinetic term that needs convincing physical explanation. 



IV. BULK GEOMETRY 



When the bulk consists of matter, the most general bulk metric for which a cosmological (FEW) metric on the 
brane is recovered, is given by a radiative black hole which is a 5-dimensional generalization 0, 0, ^| of Vaidya 
black hole In terms of transformed (null) coordinate v = t + J dr/f, the bulk metric can be written as 

dS 2 = -f(r, v) dv 2 + 2dr dv + r 2 dY,j (4.1) 
where S3 is the 3-space. For a spatially flat brane, the function f(r, v) is given by 

/(r , .,_£_!*> (42, 

with the length scale I related to the bulk (negative) cosmological constant by A5 = — 6/Z 2 and m(v) is the variable 
mass of the Vaidya black hole. This bulk black hole metric is a solution of the 5-dimensional Einstein equation with 
the bulk energy-momentum tensor 

T^n = ^QmQn (4.3) 

where qM are the ingoing null vectors and ip — dm/ dv is the rate of incoming radial energy flow to the black hole. 
Now, the black hole mass is the rescaled Weyl parameter which is further related to the scale factor by Eq (|3.2|) 

K 2 

m(t) = —C(t) cx a a (t) (4.4) 

Hence, with the scale factor of Eq l|3.8[) . the on-brane mass of the bulk black hole turns out to be 

m{t) = m Q t 2a/ ^ w+1} (4.5) 
where mo is the black hole mass at the onset of matter-dominated era, which is given by 

2/o x2q/3(™ + 1) 

m Q = C Q -[-H \ (4.6) 

Since t is the proper time on the brane, Eq (|4.5I) gives the on-brane mass m(t). In order to obtain the black hole 
geometry from the point of interest of a braneworld observer, we have to find out the off-brane mass m(v). However, 
an exact expression for the bulk geometry can never be obtained purely from the brane data but an approximate 
expression for the same at the vicinity of the brane can be obtained by following a perturbative brane-based approach. 

At low energy, the Friedman brane, located outside the "event horizon" , moves radially in the bulk. Its radial 
trajectory is given by the geodesic r(t), which reduces to the scale factor a(t) at the brane-location. To a braneworld 
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observer, a brane moving radially outwards in the bulk is identical to an expanding brane [28j. Hence, the function 
fir, v) at the brane-location reduces to 



fir, V)|br 



Co 



r a (t)\x 



Co 



a a (t) 



(4.7) 



where r(t)|b ra ne is the radial trajectory at the brane-location. Thus, the null coordinate v turns out to be 



v = t + 



I 



In 



r/l- 



a/3(w+l) 



r/l + v /m "W 3 ( tu + 1 ) 



tan 



r/l 



(4.8) 



The off-brane mass m(u) at the vicinity of the brane can be found out by expanding m{v) in Taylor series around 
its on-brane value m{t) as 



dm 




f dr 


i 


d 2 m 






dti_ 


t 1= t J 


7" 


h 2 


[ 94 _ 


t 1= t 1 


It) 



m(v) — m(t) + 



Hence in the present scenario, m(u) can be well approximated as 



l(v) « m t 2 ^ 3 ^ + Jp^ t -l+a/3(v>+l) 
3{w + 1) 



1 (r/l- V™o* q/3( "' +1) 

2 n { r/l + ^/m^t a /^ w + 1 ) 



tan 



r/l 



(4.9) 



(4.10) 



One can now substitute the above expression for m(w) into Eq 1|4.2I1 to get the function f(r, v) near the brane. It 
is straightforward to construct the relevant bulk metric at the vicinity of the brane by using together the function 
fir, v) and the coordinate v given in Eq (|4.8p . To avoid the lengthy terms involved, we skip the final expression 
for the metric. From the physical point of view, it is not required either. The null coordinate v and m(i>) given in 
equations (|4.8|) and (|4.1U|I respectively will carry all the informations about the bulk geometry to the brane. 



V. SUMMARY AND DISCUSSIONS 



In this article, we have shown that the braneworld cosmology with bulk matter can explain structure formation. 
When the bulk is constituted of matter, then the effective Einstein equation on the brane gives rise to a quantity that 
can act as an additional perfect fluid. This so-called "Weyl fluid" is the combined effect of the bulk Weyl tensor and 
bulk energy-momentum tensor, projected onto the brane. We have shown that the nature of the Weyl fluid depends 
on the energy-exchange between the brane and the bulk so that for strong bulk-brane coupling, it can dominate over 
ordinary matter. Further, following Newtonian analysis of gravitational instability, we have shown that this Weyl fluid 
can account for the required amount of fluctuations in order to explain structure formation. Thus we conclude that 
the Weyl fluid can mimic dark matter in structure formation with the advantage that nowhere we need to introduce 
any ad hoc extra matter such as dark matter. We have also investigated the bulk geometry given by a radiative Vaidya 
black hole and obtained the geometric quantities relevant for a braneworld observer. 

Throughout this article, we stick to the perturbations from the Newtonian gravitational instability. Here in no way 
we tried to study in details the braneworld perturbations as done for matterfree bulk scenario in numerous papers 
HH l30l I3TL l32| . Rather, we tried to see if the Weyl fluid can mimic dark matter in cosmological context with the 
simplest yet logical analysis of gravitational perturbations. The detailed study of scalar, metric, curvature, vector and 
tensor perturbations as well as theoretical studies of CMB anisotropics with the Weyl fluid behaving as dark matter 
are left for future works. 

An interesting issue is to study gravitational lensing that serves as a probe of structures. Significant difference 
in the bending angle due to the difference in the potentials has been reported in [20I l2l| . It is to be seen how the 
effective potential affects cosmological lensing. A comparative study of the lensing effects with those of the other dark 
matter models, e.g. scalar field haloes |33| and verification from |34j | will reveal which properties of dark matter can 
be reflected by Weyl fluid. 

Finally, we have analyzed structure formation for a flat universe without considering any accelerated scale factor. 
Having realized that braneworld gravity can account for ~ 30% of cosmic density usually attributed to dark matter, 
one can pay attention to another ~ 70% which is considered to be dark energy with negative pressure. There currently 
exists some braneworld models of dark energy [35| consistent with supernova data. We expect that this perturbative 
analysis can be applied to those models, thereby providing an unified description of dark matter and dark energy 
from braneworld gravity, with excitingly new features. 
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